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I. ABSTRACT 

Ferroic transition metal oxides, which exhibit sponta- 
neous elastic, electrical, magnetic or toroidal order, ex- 
hibit functional properties that find use in ultrastable 
solid-state memories to sensors and medical imaging 
technologies. To realize multifunctional behavior, where 
one order parameter can be coupled to the conjugate field 
of another order parameter, however, requires a common 
microscopic origin for the long-range order. Here, we for- 
mulate a complete theory for a novel form of ferroelectric- 
ity, whereby a spontaneous and switchable polarization 
emerges from the destruction of an antiferroelectric state 
due to octahedral rotations and ordered cation sublat- 
tices. We then construct a materials design framework 
based on crystal-chemistry descriptors rooted in group 
theory, which enables the facile design of artificial oxides 
with large electric polarizations, P, simultaneous with 
small energetic switching barriers between +P and -P. 
We validate the theory with first principles density func- 
tional calculations on more than 16 perovskite-structured 
oxides, illustrating it could be operative in any materi- 
als classes exhibiting two- or three-dimensional corner- 
connected octahedral frameworks. We show the princi- 
ples governing materials selection of the "layered" sys- 
tems originate in the lattice dynamics of the A cation 
displacements stabilized by the pervasive BC>6 rotations 
of single phase ABO3 materials, whereby the latter dis- 
tortions govern the optical band gaps, magnetic order 
and critical transition temperatures. Our approach pro- 
vides the elusive route to the ultimate multifunctionality 
property control by an external electric field. 



II. INTRODUCTION 

In the search for new classes of multifunctional ma- 
terials, the design or discovery of ferroelectrics in which 
the spontaneous electrical polarization couples strongly 
to other structural, magnetic, orbital, and electronic de- 
grees of freedom is a challenge being actively pursued as a 
means to achieve electric field-controllable emergent phe- 
nomena such as ferromagnetism 1,2 . Much of the current 
materials-by-design effort has focused on the structurally 
and chemically complex ABO3 perovskites, a large class 
of functional materials that display a wide range of prop- 
erties due to their highly tunable ground states. Because 
of the high susceptibility of perovskite materials towards 



polar structural instabilities, a notion has emerged that 
it is generally more productive to start with a material 
that displays, for example, ferromagnetism, and devise a 
way to induce ferroelectricity. Two highly successful ap- 
proaches that have captivated the attention of researchers 
over the last decade are that of epitaxial strain engineer- 
ing (strain-induced ferroelectricity, mutliferroicity) and 
that of selective chemical substitution of a stereochemi- 
cally inactive cation with a lone-pair-active cation such 
as Bi 3+ , as in BiFeOs 3,4 . While highly successful at cre- 
ating new multiferroics (materials that are both ferro- 
magnetic and ferroelectric), generally speaking these ap- 
proaches have not led to a widespread solution to the 
central problem of strong coupling between the polariza- 
tion and the magnetism. The reason is believed to be 
due to the common nature of the ferroelectricity in these 
materials (small cation displacements such as those in 
the prototypical perovskite ferroelectric BaTiOs). New 
ideas to realize ferroelectricity are clearly needed. 

Although perovskites are often what comes to mind 
when discussing oxide ferroelectricity, the overwhelming 
majority of oxide perovskites - particularly those which 
have active electronic, magnetic, and orbital degrees of 
freedom - adopt highly distorted, non-polar, ground state 
structures in which the B06 octahedra are rotated about 
one or more of the crystal axes 5 . A fascinating question 
that has only recently been considered in earnest, start- 
ing with the work of Ref. 6, concerns how to directly 
control these octahedral rotations with an external elec- 
tric field. Since rotations of the B06 octahedra couple 
strongly to other properties 7 10 , magnetism 11 for exam- 
ple, gaining control over them could solve the strong cou- 
pling problem 12 . Another way of stating the challenge is, 
"how can octahedral rotations induce a spontaneous po- 
larization?" 

By themselves octahedral rotations cannot induce fer- 
roelectricity in simple perovskites, but recent work has 
demonstrated that they can induce ferroelectricity in 
layered A-site ordered double perovskites and superlat- 
ties 6 ' 13 16 and in Ruddlesden-Popper phases 17,18 (and 
other layered perovskites 19 21 ). The most common real- 
ization of this novel rotation-centric, ferroelectric mech- 
anism has been referred to as hybrid improper ferro- 
electricity (HIF) 18 . The defining feature of HIF is a 
symmetry-allowed trilinear coupling in the free energy, 



^tri = 1QR1QR2P, 



(1) 



where P is the amplitude of the polarization, Qm and 
Qr2 are the amplitudes of non-polar, symmetry inequiv- 
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FIG. 1. (a) The paralectric Pnma structure (Glazer rotation 
pattern a~a~c + ) of the ABO3 constituents, (b) the ferroelec- 
tric (A/ A )B2 06 structure with the same a~a~c + rotation 
pattern, (c) the a _ a _ c° and (d) the a°a°c + rotation patterns 
that make up Pnma. (e) First principles calculated octahe- 
dral rotations of a suite of materials that span a wide range 
of tolerance factor. 

alent octahedral rotation modes 22 , and 7 is the coupling 
coefficient. The presence of this invariant in the free en- 
ergy implies that when Qr\ and Qr2 become non-zero, 
a polarization will be induced in the ground state struc- 
ture, even in the absence of prototypical, BaTi03-like, 
polar (zone-center) instabilities. 

III. SIMPLE HEURISTIC DESIGN RULES 
FROM FIRST-PRINCIPLES DATA 

A. Recent work on rotation driven ferroelectrics 
by design and a problem 

Using a combination of first-principles electronic struc- 
ture methods and symmetry arguments, Rondinelli and 
Fennie recently established design rules for the creation 
of hybrid improper ferroelectrics in ABOs/A'BOs super- 
lattices 15 . The aforementioned superlattices correspond 



to (A/A / )B2 06 double perovskites in which the A and A' 
sites order into alternating layers along [001] (note we will 
use these two different views of the structure interchange- 
ably). These rules can be summarized as follows: (1) the 
chemical criterion states that ABOs/A'BOs superlattices 
allows for a trilinear invariant, .Ftrij in the free energy by 
symmetry, and (2) the energetic criterion states that at 
least one of the perovskite constituents of the superlattice 
should have a strong tendency towards the Pnma per- 
ovskite structure (that is, Pnma should be the ground 
state structure, preferably, or a met ast able phase with a 
wide stability window). The majority of perovskites form 
in the Pnma space group 5 ' 23 ' 24 and therefore the rules 
of Rondinelli and Fennie are widely accessible to many 
chemistries and have the potential to lead to new classes 
of multifunctional materials. Note that the symmetry 
of the Pnma structure is established by two symmetry- 
inequivalent octahedral rotations, which in Glazer nota- 
tion are: a _ a _ c° (with amplitude Qr because this ro- 
tation pattern transforms like the irreducible representa- 
tion (irrep) of Pm3m) and a°a°c + (with amplitude 
Qm-, which transforms like M3"), shown in Figures 1 (c) 
and (d) respectively. The energetic criterion ensures that 
the combined rotation pattern a~a~c + survives in the or- 
dered double perovskite. 

These design rules establish when such a hybrid im- 
proper state should exist. The microscopic mechanism 
responsible for the polarization, however, is still unknown 
and hence the design rules do not directly address the 
question of whether or not this hybrid improper state is 
a functional ferroelectric or simply a pyroelectric. With- 
out insight from the microscopies, a fundamental materi- 
als problem that prevents widespread realization of these 
new multifunctional materials remains - that of under- 
standing how to design a material with a large sponta- 
neous polarization and a low ferroelectric switching bar- 
rier. This problem is best described by considering: 

Conjecture 1. Given the form of the trilinear 
coupling term, Eq. 1, it is reasonable to assume 
that the spontaneous polarization P will increase 
as the strength of each rotation, Qr\ and Qr2, in- 
creases. That is, as the energetic criterion becomes 
increasingly satisfied. Indeed, to lowest order the 
polarization can be shown to be proportional to 



where Ap is the polar mode stiffness renormalized 
by rotations and other structural distortions (see 
Appendix for details). 

Corollary 1. If this assumption is true there is 
a problem, in that large rotations necessarily lead 
to a large ferroelectric switching barrier (barring 
pathologically flat energy surfaces), since in order 
to switch P from say up to down, you must switch 
the sense of one of the BC>6 octahedral rotations. 
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Are these two points really true, and if so, is there a way 
around them in order to realize a high P, low switching 
barrier rotation-driven ferroelectric? 



B. Materials Suite to Test Conjecture 

One way to control the magnitude of octahedral ro- 
tations in perovskites is by chemical substitution of A 
or B-site cations with different ionic radii, ta and/or rp- 
This effectively allows one to control the tolerance factor, 

r = —= , (3) 

v2(t\b + r ) 

a geometric descriptor that correlates with the stability 
of a particular ABO3 material towards octahedral rota- 
tions 5 . In general, as r decreases, the greater the suscep- 
tibility of a material towards octahedral rotations. We 
consider three families of perovskites and within each 
family we consider A = Mg, Ca, Sr, and Ba, which as 
shown in Figure 1(e), span a wide range of tolerance 
factors and therefore have a wide favorability towards 
Puma. All exist in nature in the perovskite structure 
except MgZr0 3 and MgSn0 3 . 

In Figure 2(a) we plot a measure of the expected mag- 
nitude of the a~a~c + rotation pattern in the A/ A' double 
perovskites versus average tolerance factor, r avg (defined 
as the average tolerance factor of ABO3 and A'BOs). 
This measure, based on the properties of the Prima per- 
ovskites, is the geometric mean of Qm and Qr, where 
both quantities are defined as, 

0M^((Qr 3 )+(of° 3 >) 5 (4) 

where (Q^f° 3 ) and (Qm B ° 3 ) are the amplitudes of M 3 + 
in Prima ABO3 and A'BOs respectively. A nearly linear 
behavior can be seen, and as expected, the smaller the 
average tolerance factor, the larger the magnitude of the 
rotations. 5 

In Figure 2(b) we plot the energy difference, AEp, be- 
tween the lowest energy paraelectric state (the fully re- 
laxed a _ a _ c° structure in space group Pmraa), and the 
polar state in the actual ABO3/ A'BOs system. By us- 
ing this paraelectric state as a reference structure this 
not only gives us an idea as to the 'stability' of the hy- 
brid improper ground state but also the intrinsic ener- 
getic barrier to switch the polarization - via switching the 
a°a°c + rotation - between ferroelectric domains. Across 
our entire test suite of materials of sixteen different su- 
perlattices (SLs), we find that the stability of the hybrid 
improper state smoothly increases as the average toler- 
ance factor decreases, consistent with the energetic cri- 
terion of Rondinelli and Fennie, which subsequently im- 
plies that the intrinsic barrier to switch the polarization 
also increases. Therefore, in order to have a ferroelectric 
that is switchable in an experimentally realizable elec- 
tric field the average tolerance factor of the perovskite 



constituents making up the superlattice will need to be 
large. Does this mean that the spontaneous polarization 
will necessarily be small, as alluded to by Corollary 1? 

We now consider how the induced polarization varies 
with T avg . In Figure 2(c) we plot the calculated polariza- 
tion of the various superlattices versus average tolerance 
factor where it is seen that no obvious trend exists across 
all sixteen SLs. In fact, it appears to be maximal near 
~0.98. The lack of a clear correlation between P and the 
magnitude of the rotations is in stark contrast to the ex- 
pectation based on the trilinear coupling picture. Indeed 
from Equation 2 and our first-principle results, which in- 
dicate that QmQr oc (1 — T avg ) 2 as shown in Figure 2(a) 
and that Ap is largely independent of tolerance factor (as 
expected for the renormalized value in Prima as shown 
in Fig. Sll), the trilinear coupling picture would predict 

D QmQr / 1 n2 / k \ 

Ptri OC = OC (1 - T avg ) , (5) 

A P 

where the renormalized force constant of the polarization, 
Ap, is independent of r. 25 Figure 2(c), however, makes 
clear that there is not this (or any) simple, universal, 
chemistry-independent correlation between the polariza- 
tion and the magnitude of the rotations: 

P 96 (1 - r avg ) 2 . (6) 

What went wrong, if anything, in this argument? Per- 
haps the expected magnitude of the rotations in the su- 
perlattices is not approximated well by the averages? 
While this may in fact lead to a small discrepancy, it 
is unlikely that this is the origin of the qualitative dis- 
agreement between the expected polarization within the 
trilinear coupling picture and one's own chemical and 
physical intuition. Instead, the key is to realize that 
the coupling coefficient, 7, can have a nontrivial, non- 
monotonic dependence on the particular ABO3 materials 
making up the superlattice. In order to regain an intu- 
itive and predictive relationship between the polarization 
and simple chemical descriptors such as a tolerance fac- 
tor, however, the microscopic origin of the coupling con- 
stant needs to be addressed. Does 7 really depend on 
the specific material constituents of the superlattice or 
can it be described by a universal functional form that is 
independent of chemistry? 

To address this question, let us first look more closely 
at the polarization within a specific family of compounds, 
the germanate, AGeOa/A'GeOs SLs. As shown in Fig- 
ure 2(c), the polarization initially increases as r de- 
creases, but then remains fairly constant as Tavg 

further 

reduced, but then increases substantially before steadily 
decreasing as average tolerance factor is still further re- 
duced. As already concluded, this is not the r avg depen- 
dence expected from Equations 5 and 6. 

If, however, a similar exercise is performed but this 
time keeping the chemistry of one of the A-sites fixed, a 
pattern begins to emerge. This is best seen by consider- 
ing two cases. First, consider the AGeOs/BaGeOs SLs 
as the average tolerance factor decreases: 
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FIG. 2. ABOs/A'BOs superlattices: (a) Expected strength of a~a~c + rotations based on properties of ABO3 constituents, 
(b) Stabilization energy (Note: in a proper ferroelectric transition with an order parameter Q, AE oc Q. Since in our model 
Q ^ r fv g we plot AE as linear in rf vg ), (c) Polarization versus average tolerance factor with stannates highlighted, and (d) 
Polarization versus the tolerance factor renormalized by the tendency for A-site displacements showing an almost perfect fit to 
the first principles data. 



(1) SrGe0 3 /BaGe0 3 , r avg = 1.07, P=0/iC/cm 2 

\ 

CaGe0 3 /BaGe0 3 , r avg = 1.04, P=4/iC/cm 2 

\ 

MgGe0 3 /BaGe0 3 , r avg = 0.99, P=16/iC/cm 2 . 

Here the polarization monotonically increases, consistent 
with one's intuition. Next consider the MgGeOs/A'GeOs 
SLs: 

(2) MgGe0 3 /BaGe0 3 , r avg = 0.99, P = 16/iC/cm 2 

\ 

MgGe0 3 /SrGe0 3 , r avg = 0.96, P = 7/iC/cm 2 

\ 

MgGe0 3 /CaGe0 3 , r avg = 0.93, P = 2/iC/cm 2 

where even though the average tolerance factor is still de- 
creasing, in contrast to case (1) the polarization steadily 
decreases. What do we learn from this? 

Notice that in the first case the largest alkaline-earth 
cation we considered (Barium) is keep constant while the 
A-site cation steadily become smaller. In the second case, 
however, the smallest alkaline-earth cation we considered 



(Magnesium), is keep constant while the size of the A- 
cation is again steadily decreased. It would appear, at 
least within a family of compounds (for which the B- 
cation is keep constant), that the polarization is domi- 
nated by the difference in the ionic radii of the A and 
A'-cations. By elucidating the microscopic mechanism 
that leads to the trilinear coupling we will indeed argue 
that 

7 oc At U - -r^—) (J) 
\ 1 - r avg J 

where At is the difference in tolerance factors of the 
AB0 3 constituents 26 ' 27 , while 71 and 72 are coupling co- 
efficients originating from the fact that there are actually 
two distinct contributions to the polarization. This leads 
to a simple, chemically intuitive, universal result for the 
polarization: 

P a AtT [(1 - r avg )] « Ar(l - r avg ) (8) 

where J 7 is a smooth function of (1 — r avg ) and is ap- 
proximately linear for all systems except those with very 



5 



(a) 




(c) 



^ 6 r 



< 



Increasing tolerance factor (r) 



I Qr, a a c° 



I I 



o o o 

?H d CU £_. r 1 i_i rH 

|s s -a I «§ I 

s ^ ^ o o co 



O O O O Q O 



o pq 



CO CO CO 

o o o 

fl CD 03 

Q o 

PQ co pq 



FIG. 3. (a) mode, (b) for each ABO3 constituent: the stiffness of the antipolar Q x + mode of Pm3m. The ABO3 cubic 
perovskite structure is unstable to the distortion when A < 0. (c) Amplitudes of the modes in the first-principles Pnma structure 
of ABO3 materials arranged by tolerance factor. In the Pnma structure the octahedral rotations induce A-site displacements 



that transform like the irrep in Pm3m. 



small tolerance factors (so small that either the per- 
ovskite phase doesn't form or the barrier to switching 
is so large that they are no longer good candidates for 
functional materials) for which the quadratic term must 
be considered. 

Now a clear correlation between the expected polariza- 
tion and simple chemical descriptors - T avg and At - is 
realized as shown in Figure 2(d) where a nearly perfect 
fit over the entire tolerance factor range is seen. This 
result implies that it is possible to simultaneously de- 
crease the ferroelectric switching barrier, by increasing 
r avg , and increase the spontaneous polarization by in- 
creasing the difference in tolerance factors. Why? To 
answer this question we need to understand the origin of 
this rotation-induced ferroelectricity. 



IV. THEORY OF THE PROPOSED DESIGN 
RULES 

A. Rotation-driven antiferroelectricity in Pnma 
perovskite 

Turning our attention back to the ABO3 perovskite 
for a moment, it is well-known that rotations in per- 
ovskites are driven by the coordination preferences of the 
A-site cation 23 ' 24 . The combination of octahedral rota- 
tion modes that establish the symmetry of the Pnma 
structure - and R± - also favor displacements of the 
A-site cations. It turns out, because of the three dimen- 
sional connectivity of the perovskite lattice, the A-sites 
displace in an antiferroelectric pattern, the exact motion 
of which can be thought of as local, polar displacements 
confined to the two-dimensional AO layers but arranged 
180° out of phase along z, as shown in Figure 3(a). We 
denote the amplitude of this distortion as Q x + since it 

transforms like the irrep X$ in Pm3m. When we decom- 
pose the fully relaxed Pnma structures into symmetry- 
adapted basis functions of Pm3m, we find that Q x + is 
non-zero only when both Qm and Qr are non-zero, as 



shown in Figure 3(c). The mode itself, however, 
need not be unstable to appear in the Pnma structure. 
Indeed, except for compounds with very small (for per- 
ovskites) tolerance factors (r < 0.91), the test suite of 
materials considered do not favor this motion in the ab- 
sence of the a~a~c + type of rotations, as indicated by 
the positive value of the force constant, A x +; see Fig- 
ure 3(b). The appearance of a finite Q x + in the Pnma 
structure can be accounted for phenomenologically by a 
trilinear coupling in the free energy of Pm3m, 



MRXJ 



PQmQrQ j 



(9) 



This implies that once the rotations become non-zero, a 
finite antiferroelectric structural distortion, is in- 



duced 



28 



r+ (X 



QmQr 



(10) 



The similarity between the trilinear coupling of rotations 
and the antiferroelectric distortion in the ABO3 per- 
ovskites and the trilinear coupling of rotations and the 
polarization the ABOs/A'B03 perovskite SLs is curious. 
Could hybrid improper ferroelectricity actually originate 
from antiferroelectricity? 

In Figure 4(a) we plot a linear approximation to the 
layer-resolved polarization calculated from first principles 
(see Ref. 29 for an exact method) for Pnma SrSnOs 30 
where the antiferroelectric order can clearly be seen. 
With this picture in mind it is now clear that replac- 
ing alternating AO layers with A'O - for example, re- 
placing alternating SrO layers in SrSnOs with BaO to 
form (Sr/Ba)Sn206 - creates chemically inequivalent A- 
sites. The small noncancellation of the layered polar- 
ization induced by the nominally antipolar displace- 
ments results in a macroscopic polarization and a polar 
space group (Pmc2i), as shown in Figure 4b. This non- 
cancellation is the origin of HIF in this class of materials 
and a route to turn the vast number of perovskite Pnma 
antiferroelectric materials into functional ferri- electrics. 
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(a) • • • 




FIG. 4. (a) Pnma SrSnC^ and the contributions to the polar- 
ization, P = ^Player = 0, by symmetry, (b) (Sr/Ba)Sn206 
superlattice where cancellation of Player is not exact. 



We now formulate a quantitative theory of HIF that 
relies only on the properties of bulk ABO3 and A'BOs, 
and use it to derive the rule proposed above and displayed 
in Figure 2(d) that allows for the simultaneous design of 
a large P and low FE switching barrier. 



B. Cation-ordering and the design of high P, low 
switching HIFs. 

Notice above that we implicitly discussed the origin 
of the spontaneous polarization in this class of materi- 
als in the language of symmetry-lowering distortions of 
the cubic, 5-atom perovskite structure, that is, A-site 
cation ordering, antiferroelectric displacements, and of 
course the a~a~c + rotations. A natural reference struc- 
ture about which to model the polarization in the SLs is 
therefore the cubic perovskite. This has been the essen- 
tial problem with previous discussions of HIF: all have 
used the ten atom PA/mmm cation-ordered perovskite 
as the reference structure (the "problem" will become 
clear shortly). 

Imagine instead that one was to write down an expan- 
sion of the total energy in terms of the symmetry-adapted 
modes that describe the transition from the cubic Pm3m 
structure to the orthorhombic Pnma structure with the 
additional stipulation that the A-site is dynamically oc- 
cupied by two atoms, A and A 7 , with equal probability 
(hence we are considering the reference structure as a dy- 
namically disordered phase. To be able to describe the 
layered SL we therefore have to introduce an order pa- 



rameter to account for the spontaneous ordering of A and 
A' cations into layers. Group-theoretic analysis shows 
that this order-disorder transition, shown in Figure 5(a), 
lowers the symmetry of Pm3m to PA/mrnm through a 
'composition mode' transforming as X% . Further appli- 
cation of group-theoretic techniques shows that in addi- 
tion to the invariants that describe the Pm3m — > Pnma 
transition (all the symmetry allowed couplings of Qm, 
Qr, and Q x +), two nontrivial invariants are introduced 
into the free energy of Pm3m: a quadrilinear term, 

^mrx-p = JiQmQrQ x -P, ( n ) 

(where the tildes refer to modes and coupling constants 
of the 5-atom, disordered Pm3m structure) and an ad- 
ditional trilinear term 

^xix-p=l2Q xt Q x -P. (12) 

We emphasize two points. First, if not for a finite equilib- 
rium value of Q x - = (Q x -), these nontrivial couplings 
between Pnma distortions and the polarization would 
be symmetry forbidden and therefore hybrid improper 
ferroelectricity would not possible. We show this explic- 
itly from first-principles calculations of P in the a~a~c + 
structure as a function of cation-ordering, within 
the virtual crystal approximation. This first-principles 
result, shown in Fig. 5(b), shows that as the A-site or- 
dering fully saturates the polarization is maximized, val- 
idating our interpretation derived from symmetry argu- 
ments. Second, our analysis reveals that there are in fact 
two contributions to the total polarization: the first 

Pi oc 71 (Q x -)QmQr: (13) 

originates from the coupling to rotations (Eq. 11) and is 
the contribution to P usually thought of when discussing 
hybrid improper ferroelectricity. The second 

P2 0CT2(Q X -)Qx+, (14) 

originates from a direct coupling to the antipolar dis- 
placements, Eq. 12 (a third minor point is to notice that 
when X3" condenses, e.g., when one uses the 10-atom 
ordered system as a reference structure, the unit cell 
doubles along [001] and the 72 coupling causes the polar 
mode to mix with the antipolar mode resulting in a loss 
of distinction between ferroelectric and antiferroelectric 
distortions). 

Rather than just the direct rotation contribution to 
the polarization, a surprising consequence of our analysis 
shows that in general there is a second independent con- 
tribution originating from antiferroelectricity, suggesting 
hybrid improper ferroelectricity could exist without rota- 
tions. The identification of a material with rotation-less 
hybrid improper ferroelectricity would be an exciting new 
avenue to pursue multifunctional materials as the switch- 
ing barriers would be expected to be much lower than in 
proper ferroelectrics. 
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FIG. 5. (a) Layered A-site cation ordering mode, (b) Berry 
phase polarization as a function of cation ordering calculated 
within the virtual crystal approximation for (Sr/Ba)Sn2 06- 
Note Qx3- — =^ A-site is fully disordered (Pnma). 



Antiferroelectricity without rotations, however, is rare 
(if it exists at all) in perovskites. Indeed, as we discussed 
in a previous section, antiferroelectricity in Pnma per- 
ovskites is induced by the a~a~c + rotations. In such 
classes of materials our analysis may therefore be rewrit- 
ten in the form of an effective trilinear coupling. We 
integrate out the antiferroelectric mode 



ix 5 



QmQr 



72 A 



A x 5 



(15) 



to lowest order (where A Xb is the force constant of the an- 
tiferroelectric mode renormalized by rotations and other 
distortions), and minimize over Q x - to obtain 

^td = iQmQrP (16) 



with 7 ( -727^- 4 
tails), which leads to 



71) (Q x -) ( see Appendix for de- 



7" 



(Qm)(Qr) 



(17) 



perovskites we have been discussing), the conjecture that 
the polarization is simply proportional to the amplitude 
of the rotations is incorrect. Although interesting, the 
consequences of this are not profound as any physically 
reasonable force constant would be a smooth function of 
the tolerance factor and therefore cannot be the origin of 
the discrepancy displayed in Figure 2(d). 

What is interesting, is the only remaining unknown 
- the cation ordering (Q x ~). One normally thinks of 
this as a number between and 1, reflecting the crys- 
tallographic difference in occupancy of a site. Here this 
occupancy difference is irrelevant in the sense that in the 
ordered state of the Pnma SLs, the important detail (in 
relation to P) is the difference in the 'susceptibility' of 
the A and A' cations to displace from their ideal Pm3m 
positions 31 . This can be made clear by considering the 
following example. Imagine two different A/ A' cation- 
ordered materials that are both 100% ordered, but in 
one case the tendency of the A and A' cations to off- 
center is similar, whereas in the second system the A 
and A' tendency to off-center varies drastically. Even 
though cation ordering is 100% in each of these systems, 
the consequences of this ordering should differ consider- 
ing the ferrielectric mechanism described. As far as the 
polarization is concerned it is irrelevant that A and A f 
may be different atoms only to have their tendency to off- 
center be similar as this would lead to a near complete 
cancelation of the layered dipoles (just as if they were the 
same atom). In this case, it is therefore fruitful to think 
of the difference in the ability to displace between two 
A-sites as an absolute magnitude of the cation ordering. 

We propose that the absolute magnitude of the cation 
ordering is reflected in the difference in magnitude of the 
antiferroelectric mode, observed in the ABO3 per- 

ovskite from that observed in A'BOs, 

(Q x -) oc ({Q x f>) - {Q A X ?° 3 )) = A Qxs+ . (18) 

Since Q x + 1S nnear in 1 — r (see appendix), 

A Qx5+ ex At (19) 

becomes a simple crystal chemistry rule to understand 
how the cation ordering contribution to P is related to 
the difference in tolerance factor, At, between ABO3 and 

A / B Q 3 26,27 

Assuming that the QnS in the disordered perovskite 
can be approximated by the average of the respective 
QnS in each constituent, i.e., (Q n ) ~ Q n , 



This analysis makes clear that the trilinear coupling 
coefficient, 7, is composed of two factors: the magnitude 
of the cation ordering, (Q x -) and a term that goes like 

— l/(Ax 5 — 1)- Notice that the latter term differs from 
1 due to the direct coupling of the antipolar and polar 
distortions. Therefore, if A Xb has a dependence on toler- 
ance factor, then even in systems where antiferroelectric- 
ity by itself is stable (such as the majority of the Pnma 



7 oc (-72 



l-r ; 



■71) At 



(20) 



and (Qm)(Qr) oc (1 — r avg ) 2 we therefore have 



P = -^i{Qm)(Qr) 
Ap 



oc (-72 



1 " T K 



- 7 i)Ar(l -r avg ) 2 



(21) 



8 



ciQmQr 

C2Q X +* 




-60 



0.85 



0.90 



0.95 



1.00 



1.05 



1.10 



T 3 



avg 



(b) 



600 



300 



-300 



-600 




0.85 



0.90 



0.95 



1.00 



1.05 



1.10 



avg 



FIG. 6. (a) Polarization contributions from fitting to first principles mode amplitudes scaled by the difference in magnitude 
of the antiferroelectric mode, (b) polarization contributions from tolerance factor model scaled by the difference in tolerance 
factor. 



and therefore 

P = At (ci(1 - r avg ) 2 + c 2 (l - r avg )) (22) 

which shows a nearly perfect fit to the first-principles cal- 
culations across all three families of compounds as dis- 
cussed previously and shown in Figure 2(c). In Figure 6 
we show each contribution to the polarization separately 
across our material suite. Each contribution, computed 
using only the first-principles structural distortion ampli- 
tudes, is a smooth monotonic function of tolerance factor 
when scaled by the difference in magnitude of the antifer- 
roelectric mode. Furthermore, the rotation-driven term 
from Equation 13 has an oppositely oriented polarization 
than the antiferroelectric-driven term from Equation 14. 
Notice that the rotation-driven polarization term dom- 
inates for small tolerance factor materials. Quite sur- 
prisingly, however, it is the antiferroelectric contribution 
that dominates for materials with r avg close to 1 where 
the switching barrier is small. 

This quantitative model results in a simple, chemically 
intuitive design rule of thumb. 

To keep the switching barrier low and the polariza- 
tion high: two bulk ABO3 and A f BOs perovskites 
should be selected such that the average tolerance 
factor, r avg , is maximized and the difference in 
their tolerance factors, At, is also large, 

and for most reasonable tolerance factors the model is 
equivalent to 

P~ Ar(l-r aV g). (23) 

C. Extension to Perovskite-derived cation ordered 
Ruddlesden-Popper materials 

The chemical intuition learned in the previous sections 
is applicable to more than just perovskite heterostruc- 
tures having a 3D network of corner-connected octahedra 



such as the 1-1 SL. For example, it has been known for 
some time that certain n = 2 Ruddlesden-Popper (RP) 
compounds - layered perovskites with disconnected oc- 
tahedra along [001] - display rotations of their octahedra 
and are polar 17 , as shown in Figure 7(a). Recently, one 
such system, Ca3Mn207 18,32 , has been shown to display 
a strong coupling between the rotation-induced ferroelec- 
tricity and magnetism. An analysis by Aleksandrov and 
colleagues has suggested that the a~a~c + octahedral ro- 
tation pattern is responsible for the polar state. Is this 
really true? 

Using our suite of ABO3 perovskites we studied nine 
A3B2O7 systems spanning a wide range of ABO3 toler- 
ance factors. In Figure 7(b) we plot the amplitudes of 
the a~a~c° and a°a°c + distortion patterns. For the six 
A3B2O7 derived from Pnma perovskites, we find that 
the RP structure displays a similar rotation pattern and 
is polar. This indeed suggests that any Pnma perovskite 
would be polar if it can be synthesized in the n = 2, 
A3B2O7 RP structure. 

In order to realize a functional ferroelectric, however, 
the polarization should be large with a low switching bar- 
rier. Based on what we learned in the previous section 
we would expect that as the tolerance factor of the ABO3 
perovskite making up the RP material decreases, the oc- 
tahedral rotations, and consequently the polarization and 
the switching barrier, would get larger. As shown in Fig- 
ures 7(b), Figure 8(a), and Figure 8(b), this is exactly the 
case. Perhaps the large barrier to switch the polarization 
in Ca3Mn207 or Ca3Ti207 (two known polar RPs) is the 
reason why ferroelectricity has yet to be shown experi- 
mentally. (Note, Sr3Sn207, which has been synthesized 
previously, 30 is an ideal system to explore for ferroelec- 
tricity). 

Is the origin of the polar state in these Pnma RPs 
similar to what we have been discussing? First note 
that unlike the perovskite superlattices, the polarization 
is a smooth and monotonic function of tolerance factor. 
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This behavior is in fact consistent with the mechanism 
we outlined in the previous sections and can be easily 
understood by examining the RP structure. First, in 
even layered RPs, such as the n = 2, the interface be- 
tween the rock salt layer and the perovskite layer breaks 
the local inversion symmetry of the individual octahedra, 
therefore cation-ordering in not required as in the per- 
ovskite super lattices. Additionally, as shown in the inset 
to Figure 8(b), the a~a~c + octahedral rotation pattern 
still induces antipolar A-cation displacements in the RP. 
Given that the AO layers at the interface and the layer 
between the perovskites are different by symmetry, a net 
polarization will arise - hybrid improper ferroelectricity 
without cation ordering. Unlike the superlattice example 
however, the noncancellation of the A-cation displace- 
ments from this effect is very small, yet there is still 
a substantial polarization because of the additional AO 
layer present in the unit cell of the RP. The polariza- 
tion varies smoothly because both of these contributions 
depend only on one type of A-cation antiferroelectric dis- 
placement and therefore only on one ABO3 tolerance fac- 
tor. 

Having established a common origin, can we use the 
chemical mismatch strategy, At, learned in the previ- 
ous section, to reduce the switching barrier (by increas- 
ing r avg ) while increasing the polarization in the A3B2O7 
structure? The design rules tell us to create A/ A' cation 
ordered RP structures, as shown in Figure 9. Tables I 
and II show the first-principles polarization and switch- 
ing barrier in the (ASn03) 2 A'0 and (AZr0 3 ) 2 A'0 mate- 
rials, respectively. Taking (CaSn03) 2 CaO as an example, 
and reading down the first column of Table I, we find 

(CaSn0 3 ) 2 CaO, P=12 /iC/cm 2 , AE=406 meV 

\ 

(CaSn0 3 ) 2 SrO, P=14 /iC/cm 2 , AE=234 meV 

\ 

(CaSn0 3 ) 2 BaO, P=17 /iC/cm 2 , AE=122 meV 



(a)< 



(a) Ferroelectric (b) 
a~a~c + A3B2O7 



Increasing ABO3 tolerance factor 




Ferroelectric ( 



\a a c 



1 1 



000000000 

^ c<3 § c3 ci 



FIG. 7. (a) The ferroelectric a~ a~c + structure of the A3B2O7 
Ruddlesden-Popper phase, (b) First principles amplitudes of 
the two octahedral rotations and the induced polar mode for 
a suite of A3B2O7 materials, arranged by increasing tolerance 
factor of the ABO3 parent. 
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FIG. 8. (a) Stabilization energy, and (b) polarization of the 
n — 2 Ruddlesden-Popper materials arranged by the toler- 
ance factor of the ABO3 end member. Inset: the A-cation 
displacement pattern in the ferroelectric a~ a~ c + structure. 



which is remarkably consistent with our superlattice de- 
sign rules. While substitution of the A' cation (in be- 
tween the perovskite blocks) in (CaSn0 3 ) 2 CaO reduces 
the switching barrier and increases the polarization, Ta- 
ble I shows that substitution of the A cation (in the rock 
salt layer) reduces the barrier but also decreases the po- 
larization. For example, reading across the first row we 
find 

(CaSn0 3 ) 2 CaO, P=12/iC/cm 2 , AE=406 meV 

\ 

(SrSn0 3 ) 2 CaO, P=3.4/iC/cm 2 , AE=182 meV 

\ 

(BaSn0 3 ) 2 CaO, P=3.5/iC/cm 2 , AE=145 meV 



Does this contradict our chemical mismatch design strat- 
egy? 

Let us first look at (CaSn0 3 ) 2 CaO where A=A . As 
Figure 10(a) shows, because of the weak symmetry break- 
ing of the perovskite/rocksalt interface, the polarization 
originating from one of the AO layers nearly cancels that 
coming from the A O layer and therefore in this material 
the polarization is approximately that originating from 
one CaO layer. Now with the substitution of the A site 
with a different cation, as shown in Figure 10(b) the AO 
layer and the A O layer no longer cancel. In the case 
where the A cation radius is much larger, for example 
A =Ba, the polarization in this A O layer is nearly zero. 
According to the model we derived in the previous sec- 
tions, this should lead to a polarization that is approxi- 
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FIG. 9. (a) The n = 2 Ruddlesden- Popper structure (with- 
out octahedral rotations) with a single A-cation. (b) A/A' 
ordered n — 2 phase, (c) A/ A' ordered n — 1 phase. All three 
of these layered perovskites are polar in the a~ a~ c + ground 
state. 
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FIG. 10. Simple model for the polarization showing the AO 
and A O layer contributions, (a) A=A , (b) A is replaced by 
a larger cation (for instance Ba replaces Ca) (c) A is replaced 
by a larger cation. 



mately twice the contribution from one CaO layer multi- 
plied by the ratio of the average AB0 3 tolerance factors, 
i.e., P « 15/iC/cm 2 , which compares amazingly well to 
the first-principles result of 17/iC/cm 2 . Going across the 
row, the tolerance factor increases much faster than the 
mismatch increase and the polarization decreases, e.g., 
the model predicts P ~2.3/iC/cm 2 (first-principles value 
equals 3.5 /iC/cm 2 ) in (BaSn0 3 ) 2 CaO because the polar- 
ization from one CaO layer is reduced by the substantial 
increase in the average tolerance factor as shown in Fig- 
ure 10(c). 

The same general design rules can be applied to the 
n = 1 member of the Ruddlesden-Popper series. Al- 
though the a~a~c + rotation pattern does not induce fer- 
roelectricity in A 2 B04 as there are an even number of A 
cations per perovskite block, the a~a~c + ground state 
of the A/A' cation ordered (AB0 3 )AO is polar 12 , for ex- 
ample: 

(CaSn0 3 )SrO, P = 5.7 /iC/cm 2 , AE = 235 meV 

\ 

(CaSn0 3 )BaO, P = 11 /iC/cm 2 , AE = 159 meV 
which is consistent with our general design rules. 



V. CONCLUSION 

We have shown how to take the most common per- 
ovskite, the Pnma perovskite, and turn it into a func- 
tional ferroelectric. We have presented the microscopic 
theory of this new type of ferroelectricity, based on the 
interplay of steric octahedral rotations, antiferroelectric- 
ity, and cation ordering, through both group theoreti- 
cal methods and systematic first-principles calculations 



on more than 16 different ordered ABO3/A B0 3 super- 
lattices and 20 different (AB0 3 ) 2 A / and (ABOs^A'O 
Ruddlesden-Popper (RP) compounds. Using results only 
on super lattices, we have derived a design strategy ex- 
pressed only in terms of crystal-chemistry descriptors 
of the ABO3 constituents. We have shown that these 
straightforward design concepts also capture the basic 
physics of the RP systems as well and are likely general, 
universal concepts. There are no chemical restrictions on 
the type of A and/or B cations that can be used, making 
this ferroelectric mechanism compatible with the types of 
correlated electron phenomena, for example magnetism, 
usually associated strongly with the Pnma perovskites. 
This materials design framework for achieving large elec- 
trical polarizations and low switching barriers will be a 
crucial guide to the successful integration of these new 
multifunctional ferroelectric materials into next genera- 
tion devices. 



TABLE I. Ground state total polarization and switching en- 
ergy barrier of the six possible (ASnOs) 2 A'0 superlattices. 
"-" indicates compounds are stable in the paraelectric struc- 
ture without any octahedral rotation about the [001] axis. 



A'O 


(CaSn0 3 ) 2 A'0 


(SrSn0 3 ) 2 A'0 


(BaSn0 3 ) 2 A'0 


CaO 


12 /iC/cm 2 
406 meV 


3.4 /iC/cm 2 
182 meV 


3.5 /iC/cm 2 
145 meV 


SrO 


14 /iC/cm 2 
234 meV 


3.9 /iC/cm 2 
37 meV 


1.8 /iC/cm 2 
10 meV 


BaO 


17 /iC/cm 2 
122 meV 
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TABLE II. Ground state total polarization and switching en- 
ergy barrier of the six possible (AZrOs) 2 A'O superlattices. 



A'O 


(CaZr0 3 ) 2 A'0 


(SrZr0 3 ) 2 A'0 


(BaZr0 3 ) 2 A'0 




16 /iC/cm 2 
430 meV 


6.0 /iC/cm 2 
192 meV 


4.4 /xC/cm 2 
129 meV 


SrO 


17 /iC/cm 2 
293 meV 


7.2 /xC/cm 2 
67 meV 


3.4 /xC/cm 2 
19 meV 


BaO 


24 /iC/cm 2 
193 meV 


11 /xC/cm 2 
10 meV 





VI. METHOD 

First-principles calculations are performed within den- 
sity functional theory as implemented in Quantum 
Espresso 33 using ultrasoft pseudopotentials and the 
PBEsol functional 34 , which provides an improved de- 
scription of structure 34 over LDA or PBE 35 . We used 
a plane wave energy cutoff of 650 eV. We used a k-point 
grid equivalent to an 8x8x8 grid of the cubic perovskite. 
The total polarization was calculated using the Berry 
phase method. 
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VIII. APPENDIX: LANDAU THEORY 

The choice of high symmetry reference structure is a 
key step in modelling materials in which multiple distor- 
tion modes couple in nontrivial ways. Prior work on hy- 
brid improper ferroelectrics used the A/ A' cation ordered 
P4/mmm structure, rather than the cation disordered 
cubic PraSra, as the high symmetry reference structure. 
In this appendix we provide the Landau theory for each 
of these reference structures. 



A. Polarization Starting from the Cation Ordered 
P4 /mrnrn Structure 

In the 10 atom, cation ordered P4/mmm structure (see 
Figure 5a) there are three key distortion modes that we 
will treat as order parameters: Qri which describes the 



strength of a~a~c° octahedral rotations, Qr2 which de- 
scribes the strength of a°a°c + octahedral rotations, and 
the polarization. The free energy around the P4/mmm 
structure is 

i ' 1 ( 24) 

^^Y,B l3 Q\Q)-r 1 Q R1 Q m P 

ij 

where the summation index denotes the three order pa- 
rameters Qr\, Qr2, and P. All of the ferroelectric ma- 
terials we consider have Ari,Ar2 < 0. In other words, 
the P4/mmm structure is unstable with respect to the 
a~a~c + rotation pattern. We set Qri = (Qri) and 
Qr2 = (Qm) to obtain the effective free energy in the 
polar a~a~c + structure. 

T = \a p P 2 + 7 (Q R1 ) (Q R2 ) P (25) 

where Ap is the polar mode stiffness after renormal- 
ization from the biquadratic couplings with (Qri) and 
(Qr2) and any other modes which are symmetry- allowed 
in the polar structure. In Figure 11 we show how the 
renormalized polar mode is stable, Ap > 0, even when 
the high symmetry structure has a polar instability. 
From Equation 25 we therefore conclude the polariza- 
tion arises from the trilinear coupling to the octahedral 
rotations. Minimizing over P leads to the lowest order 
spontaneous polarization 

p = JQriHQr2) m ( 26 ) 

Ap 

In the main text Figure 2 we show that (Qri) (Q112) oc 
1 — Tavg- Figure 11 shows that the renormalized polar 
mode force constant is roughly independent of tol- 
erance factor, and therefore from this analysis we would 
expect P oc 1 — T avg . As discussed in the main text 
this tolerance factor dependence is not observed in our 
first principles calculations. This apparent problem is 
resolved by using the A/A' cation disordered Pm3m ref- 
erence structure. 



B. Polarization Starting from the Cation 
Disordered Pm3m Structure 

Instead of the 10 atom PJ^/mmm structure we will con- 
sider the Landau free energy expanded around the 5 atom 
cation disordered Pm3m structure. The P4/mmm space 
group is then a subgroup of Pra3ra, connected by the 
irrep X% which we physically attribute to cation order- 
ing of the A-site (see Figure 5). We will now demon- 
strate how A-site ordering leads to an effective trilinear 
coupling between a~a~c + rotations and the polarization. 
The spontaneous polarization will be given by Equation 
26 except that the P4/mmm trilinear coupling, 7, will 
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be replaced by an effective coupling which will depend 
on the magnitude of the A/ A' antiferroelectric displace- 
ments. 

In the 5-atom, cation disordered Pm3m structure there 
are five relevant mode amplitudes that we take to be the 
order parameters in our Landau theory: Qm, Qr, Q x +, 

Q x ~i and P (see main text for details). To fourth order 
the free energy around the Pm3m structure is 



i i 



(27) 



i2Q x +Qx- p + iiQmQrQ x - p 



where the summation index denotes the five order pa- 
rameters specified above. 

We first integrate out the antipolar mode Q x + to ob- 
tain an effective free energy in the remaining order pa- 
rameters. We assume the rotations sterically induce this 
antipolar motion and thus A x + > (which is true for 

r > 0.91, see Figure 3). We obtain 



A' 



■QmQi 



72 



(28) 



where A v+ is the mode stiffness after renormalization 
from the biquadratic couplings of each other mode to 



and to lowest order A, 



Finally, in the 



cation ordered 10 atom cell (see the inset of Figure 3(c) 
in the main text) we set Q x - ~ (Qx 3 ~) to obtain the 
effective free energy for the three hybrid improper order 
parameters Qm, Qr and P. By substitution of Q x + and 




Ap from Pm3mA 
Ap from Pnma • 



1.10 



avg 



FIG. 11. The average polar mode stiffness (force constant) 
of ABO3 and A'B03 in the high symmetry Pm3m structure 
and renormalized in the Pnma a~ a~ c + structure. 



Q x - in the free energy of Equation SI, we obtain 




(29) 



QmQrP 




where the summation index now denotes the two rotation 
modes and the polarization, and the new effective coef- 
ficients Am-, Ar, Ap and Bm,r are modified from the 
values in Equation SI through the coupling terms. The 
effective trilinear coupling constant between Qm, Qr and 
P is 



(30) 



C. Polar-Antipolar Mixing 

In this section we formally demonstrate how the cation 
ordering of the A-site mixes the polar and antipo- 
lar modes of the constituent ABO3 perovskite. The 
cation ordering along the [001] direction (equivalent to a 
ABOa/A'BOa superlattice) is accounted for in our Lan- 
dau theory by the 'composition' order parameter Q x - 
(see main text for details). In the 10 atom cation or- 
dered cell, the free energy of the disordered perovskite 

in Equation 27 is modified by Q x - = (^Q x -^ which 
leads to an effective free energy of the remaining order 
parameters, 



i i 

+ lj2 B aQ 2 iQ 2 j + PQmQrQ x + 



(31) 



' 72 



; x +P + 7 1 (Q x -)QmQrP- 



The bilinear coupling of Equation 31 leads to a mixed 
normal mode in the antipolar and polar distortions, Q x + 

and P (see Figure 3 in the main text). We can write the 
mixed mode as 



Q EB = cos(O)P + sm(0)Q 



(32) 



The mixed mode transforms with same irreducible rep- 
resentation as the polarization in the cation ordered 10 
atom cell. The strength of mixing is determined di- 
rectly from the bilinear coupling coefficient j 2 (Q x -^ 
and the relative stiffnesses of the antipolar and polar 
A 



modes A A = | (^A x + — Apj, through the relation 
72 (Q x - ) tan (0) = AA — J AA* (q x S . 



(33) 
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